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ABSTRACT 

Here we report on the significant role of a so far overlooked dynamical aspect, namely a 
secular resonance between the dwarf planet Ceres and other asteroids. We demonstrate that this 
type of secular resonance can be the dominant dynamical factor in certain regions of the main 
asteroid belt. 

Specifically, we performed a dynamical analysis of the asteroids belonging to the (1726) 
Hoffmeister family. To identify which dynamical mechanisms are actually at work in this part of 
the main asteroid belt, i.e. to isolate the main perturber(s), we study the evolution of this fam¬ 
ily in time. The study is accomplished using numerical integrations of test particles performed 
within different dynamical models. The obtained results reveal that the post-impact evolution of 
the Hoffmeister asteroid family is a direct consequence of the nodal secular resonance with Ceres. 

This leads us to the conclusion that similar effects must exist in other parts of the asteroid 
belt. In this respect, the obtained results shed light on an important and entirely new aspect 
of the long-term dynamics of small bodies. Ceres’ fingerprint in asteroid dynamics, expressed 
through the discovered secular resonance effect, completely changes our understanding of the way 
in which perturbations by Ceres-like objects affect the orbits of nearby bodies. 

Subject headings: celestial mechanics — minor planets, asteroids: general — minor planets, asteroids: 
individual (Ceres) 


1. Introduction 


Orbital resonances exist everywhere in the So¬ 
lar System, and play an essential role in the dy¬ 
namics of small bodies. The synergy of fast and 
slow orbital angles pr oduces a great assortment of 
resonant phenomena ( Whhmnsj^jhjjlkne^ J^ 

icIlQM 


Dermott fc MurravIlQSlL Milani fc Knezevidl994 

Nesvornv fc Morbidellil Il998 ). Over the years 


numerous methods and models have been devel¬ 
oped to interpret the complex dynamical environ¬ 
ment of the main asteroid belt. It is well known 
that this region is sculpted by a web of mean- 
motion and secular resonan ces coupled with sub- 


tie non-gravitational forces (Gladman et al.l 1997 

Farinella & Vokrouhlickv 

I999I: Bottke et al.l2006 

Minton & Malhotral 20If 

Novakovic et al. 201011 


The implications of these effects on a large number 
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of examples with unique dynamical characteristics 
have already been successfully described by exist¬ 
ing dynamical models. It is, however, still not 
possible to explain or predict all the dynamics of 
the main asteroid belt. 


There are two general types of orbital reso¬ 
nances in the Solar system. The most intuitive 
type, referred to as mean-motion resonances, oc¬ 
curs when the orbital periods of an asteroid and 
of a perturber are nearly commensurate. The sec¬ 
ond type, called secular resonance, concerns slowly 
varying angles like the longitude of perihelion or 
the longitude of the ascending node. 

Secular resonances may play a significant role 
in the long-t erm dynamical stability of a plane¬ 


tary system (LaskM 19891 : iKnezevic et al.l [l9^ 
Michel fc Froeschl^l997 ), however until now stud¬ 


ies related to the dynamics of small solar system 
objects considered only planets as important per- 
turbers. 


The role of the most massive asteroids in the as¬ 
teroid dynamics is generally assumed to be small, 
and in most cases it is neglected. Still, it is known 
that perturbations arising from the most mas¬ 
sive asteroids could be important in some sit¬ 
uations. Being located relatively close to each 
other, the most important interaction between 
the most massive and the rest of the asteroids 
intuitively occurs during their mutual close en¬ 
counters. The long-term effects of these encoun¬ 
ters have been studied by many authors, typically 


ilies (Nesvornv et al. 20021: Carniba et al.l 

2003 

Novakovic et a 

. 2 OIOI: Delisle & Laskad 

2012 

Carruba et al.l 

2 OI 3 II. Moreover, it was demon- 


are populations of asteroids in the 1/1 mean mo¬ 
tion resonances with the two most massive objects 
in the main belt, namely (1) Ceres and (4) Vesta. 

Nevertheless, the importance of massive aster¬ 
oids for secular dynamics and long-term chaotic 
diffusion is generally accepted to be negligible; 
thus, it has never been studied. In this paper, 
we show that this paradigm not only lacks justifi¬ 
cation, but it is actually incorrect. The results 
obtained here reveal that a nodal secular reso¬ 
nance with (1) Ceres, namely s — Sc, plays a key 
role in the dynamics of asteroids belonging to the 
Hoffmeister family. 


2. Methods and results 


The motivation for our work comes from the 
unusual shape of t he Hoffmeister asteroid family 
l)Milani et al.ll2014 1 when projected on the proper 
orbital semi-major axis Op versus sine of proper 
orbital inclination sin{ip) plane (Fig. [T]). The 
distribution of family members as seen in this 
plane clearly suggests different dynamical evolu¬ 
tion for the two parts of the family delimited in 
terms of semi-major axis. The part located at 
Up < 2.78 AU is dispersed, and seems to un¬ 
dergo significant evolution in inclination, contrary 
to that at Op > 2.78 AU, which looks much more 
condensed and practically shows no similar evo¬ 
lution. Our goal here is to reveal the mechanism 
responsible for the observed asymmetry. 


PLACE FIGURE 1 HERE. 


Analyzing the region around the Hoffmeister 
family we found a few potentially important dy¬ 
namical mechanisms. The family is delimited in 
terms of semi-major axis by two mean motion res¬ 
onances, the 3J-1S-1 three body resonance with 
Jupiter and Saturn at 2.752 AU, and the 5/2J 
mean-motion resonance with Jupiter at 2.82 AU. 
Moreover, the region is crossed by the zi=g — 
ge + s — sq secular resonance, with g, s, ge and sg 
being the secular frequencies of the asteroid’s and 
Saturn’s orbits. 

2.1. Numerical simulatious 

2.1.1. Dynamical model 

To identify which ones, if any, among the dif¬ 
ferent possible dynamical mechanisms are actu¬ 
ally at work here, we performed a set of numer¬ 
ical integrations. For this purpose we employed 
the ORBIT9 integrator embedded in the multi¬ 
purpose OrbFit packag^. The dynamical model 
includes the gravitational effects of the Sun and 
the four outer planets, from Jupiter to Neptune. 
It also accounts for the Yarkovsky thermal effect, 
a subtle non-gravitational force due to the recoil 
force of anisotropically emitted thermal radiation 


^Available from http://adams.dm.unipi.it/orbfit/ 
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by a rotating body (jBottke et al.l 120061) . causing 
mainly a secular drift in semi-major axis. The in¬ 
direct effect of the inner planets is accounted for 
by applying a barycentric correction to the initial 
conditions. 

Our simulations follow the long-term orbital 
evolution of test particles initially distributed ran¬ 
domly inside an ellipse determined by the Gauss 
equations. This ellipse corresponds to the dis¬ 
persion of the Hoffmeister family members im¬ 
mediately after the breakup event, assuming an 
isotropic ejection of the fragments from the par¬ 
ent body. 

The total number of particles used is 1678, 
the same as the number of asteroids we currently 
identified as members of the family. The family 
membership is determined utilizing the hierarchi¬ 
cal cluste ring method and stan dard metric as pro¬ 
posed by Zappala et al.l ( IQQOll . 

For simplicity, the Yarkovsky effect is approxi¬ 
mated in terms of a pure along-track acceleration, 
inducing on average the same semi-major axis drift 
speed da/dt as predicted from theoryH Assum¬ 
ing an isotropic distribution of spin axes in space, 
to each particle we randomly assign a value from 
the interval ±{da/dt)max, where {da/dt)max is the 
estimated maximum of the semi-major axis drift 
speed due to the Yarkovsky force. The value of 
{da/dt)max is determined us ing a model of the 
Yark ovsky effect developed bv IVokrouhlickvl ( 19981 
I 1999 II . and assuming thermal parameters appro¬ 
priate for regolith-covered C-type objects. In par¬ 
ticular, we adopt values of ps = Pb = 1300 kg in~ ^ 
for the surface and bulk densities ( CarrvI 12012 1. 


r = 250 J m ^ s K ^ for the surface ther- 


^This model of the net Yarkovsky force is a reasonable ap¬ 
proximation over short timescales, but may not be accu¬ 
rate enough in the long term, because the spin axis or 
the rotational period may change. The spin evolution of 
an asteroid not subject to collisions, is expected to be 
dominated by the Yark ovsky-O’Keefe-Ra dzievskii-Paddack 
(YORP) effect (see e.g. iRubincamI[20001) . The models pre¬ 
dict that YORP torques may evolve bodies toward asy mp- 
totic rotational states jCapek &: VokrouhlickS^ |2004|) . or 
could cause reshaping that would significantly increase the 
time over which objects can pre serve their sense of rota¬ 
tion l[Cotto-Figueroa et al.l[2015l) . This has implications 
for the Yarkovsky effect, however, the constant Yarkovsky 
drift we used here represents the long-term average of this 
effect. Regardless of the actual behavior of any single body, 
the average drift rate should be nearly constant for a large 
enough statistical sample. 


mal inertia ( Delbo fc Tanga 20^, and e = 0.95 
for the thermal emissivity parameter. In this way 
we found that for a body of H = 1 km in di¬ 
ameter {da/dt)max is about 4 x 10“^ AU/Myr. 
Next, we select sizes of the test particles equal 
to sizes of the Hoffmeister family members, esti¬ 
mated using their absolute magnitudes provided 
by AstDys databasjl, and geom etric albedo of 
Py = 0.047 ( Masiero et ahlEoilh . Finally, since 
Yarkovsky effect scales as oc l/U, the particle sizes 
are then used to calculate corresponding value of 
{da/dt)max for each particle, by scaling from the 
reference value for objects of I? = 1 km. 

The orbits of the test particles are propagated 
for 300 Myr, which is a rou gh estimate of the age 
of the Hoffmeiste r family ( Nesvornv et al 1 l200,5l: 
Sooto et'^ I 2 OI 5 II . Time series of mean orbital 


elements ^ are produced us ing on-line digital fil¬ 
tering (jCarpino et al.lllOSTI l. Then, for each par- 
ticle we compute the sy nthetic proper elements 
( Knezevic fc Milanill2000ll for consecutive intervals 
of 10 Myr. This allows us to study the evolution of 
the family in the space of proper orbital elements. 


If our first dynamical model were complete we 
should be able to reproduce the current shape of 
the family. However, we found that the shape can¬ 
not be reproduced with the afore described model. 
In particular, in these simulations we observed 
only a dispersion of the semi-major axis caused 
by the Yarkovsky effect, but no evolution in incli¬ 
nation (Fig. [5]). This implies that neither mean- 
motion nor secular resonances involving the major 
outer planets are responsible for the strange shape 
of the Hoffmeister family. 


PLACE FIGURE 2 HERE. 


In order to clarify the situation we turned our 
attention to the inner planets, and their possible 
role in the evolution of the family. For a represen¬ 
tative sample of about 200 test particles we re¬ 
peated the above described simulations using a 
model with seven planets (from Venus to Nep- 


®http: //hamilton. dm.unipi.it/astdys/ 

“^The mean orbital elements are obtained by removal of the 
short-periodic perturbations from the instantaneous oscu¬ 
lating elements. 
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tune), which also includes the Yarkovsky thermal 
force. However, the 7-planet model did not give 
any different results, which remain practically the 
same as the one obtained within the model with 
the four outer planets only. Thus, obviously there 
is still something missing in the dynamical model. 

2.1.2. Extended dynamical model 

Being left with almost no other option, we 
turn our attention to asteroid (1) Ceres. Hav¬ 
ing a proper orbital semi-major axis of 2.767 AU, 
and therefore being inside the range covered by 
the family members, it seems to be the only re¬ 
maining candidate. Hence, we again numerically 
integrated the same test particles, using the 4- 
planet dynamical model but this time also includ¬ 
ing Ceres as a perturbing bodjU. 

These new simulations already after about 
150 Myr very nearly matched the current spread¬ 
ing of the family in the semi-major axis versus 
inclination plane, clearly implicating Ceres as the 
culprit for what we see today. The striking feature 
observed in these runs is a fast dispersion of or¬ 
bital inclinations for objects with semi-major axis 
of about 2.78 AU. After 15 Myr of evolution the 
spread in sine of inclination is 3 times larger than 
the initial one, as can be seen in the middle panel 
in Fig. [2j 

2.1.3. Mechanism of Ceres perturbations 

The fact that the main perturber is Ceres raises 
a very important question. What is the exact 
mechanism by which Ceres is perturbing the mem¬ 
bers of the Hoffmeister family to such a high de¬ 
gree? The current paradigm suggests this may 
be the result of close encounters, or it might be 
the consequence of the 1/1 mean motion reso¬ 
nance with Ceres. However, our simulations un¬ 
doubtedly show that none of these two mecha¬ 
nisms could explain the evolution of the family. 
Actually, we found that most of the evolution is 
taking place within a narrow range of the semi¬ 
major axis. However, this range does not corre¬ 
spond to the location of the 1/1 resonance with 
Ceres, neither is there any reason for close en¬ 
counters to affect only objects within this specific 
range of semi-major axes. Moreover, it is very 

®In these simulations, for mass of C eres we use d a va lue of 
4.757 X 10“as estimated bv iBaer et al.l l201lf) . 


unlikely that these two effects would primarily af¬ 
fect orbital inclinations. Finally, regardless of the 
mechanism, such a large perturbation on aster¬ 
oids caused by Ceres has never been observed in 
the asteroid belt. This situation motivates us to 
test some other mechanisms which are at work in 
this region, despite generally being accepted to be 
negligible. These are the secular resonances with 
asteroid Ceres. 

Analyzing the secular frequencies of the Hoffmeis¬ 
ter family members we immediately notice that 
some of these are very close to the nodal fre¬ 
quency of Ceres (sc = —59.17 arcsec/yr). This is 
an interesting fact because the secular resonances 
involving the nodal frequencies are known to affect 
mainly the orbital inclination. 

To better understand the reasons for the disper¬ 
sion in inclination, and to determine the possible 
role of the secular resonance with Ceres, we pick a 
few particles that experienced significant changes 
in inclination during our numerical simulations, 
and analyzed their behavior in more detail. This 
analysis identified a mechanism responsible for the 
evolution of orbital inclinations, revealing a com¬ 
pletely new role of Ceres in asteroid dynamics. 

As an illustration, in Fig. [3] we show the evolu¬ 
tion of one test particle. This particle was initially 
located at a semi-major axis of 2.788 AU, with its 
Yarkovsky induced drift set to be negative, forcing 
it to move towards the Sun. After about 65 Myr 
this particle enters the region where the fast dis¬ 
persion in orbital inclination has been observed. 
During the time spent in this area, the particle’s 
inclination has experienced a fast increase, with 
the average value jumping from about 4.4 to 4.9 
degrees. 


PLACE FIGURE 3 HERE. 


Let us recall here that a resonance occurs when 
the corresponding critical angle librate^. In the 
case of the s — Sc secular resonance the critical an¬ 
gle is (T = D—He, where H and He are the longitude 
of the ascending node of an asteroid and Ceres, 


®A libration is the oscillation of an angle around a fixed 
point, contrary to the circulation when the angle cycles 
over all values from 0 to 360 degrees. 
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respectively. Clearly, as can be seen in Fig. [3l 
the period of increase in inclination exactly cor¬ 
responds to the period of libration of the critical 
angle of the vic = s — Sc secular resonance. This 
correlation is a direct proof that this resonance is 
responsible for the evolution and observed spread 
in the orbital inclination for asteroids belonging 
to the Hoffmeister asteroid family. In Fig. [2] we 
plotted the location of the vic resonance. 

Though our numerical simulations clearly show 
that passing through the vic resonance may cause 
significant changes in orbital inclination, we fur¬ 
ther investigated the mechanism. A key to under¬ 
stand the observed behavior are the cyclic oscilla¬ 
tions in inclination for objects trapped inside this 
resonance (see Fig. n. In the scenario with the 
Yarkovsky effect included in the model, some of 
the objects are reaching the border of the while 
steadily drifting in semi-major axis, and enter it at 
random value of the inclination cycle. During the 
time spent inside the resonance their inclination is 
continuously repeating the cycles. However, as the 
semi-major axis is evolving due to the Yarkovsky 
effect, the objects must sooner or later reach the 
other border of the i^ic, and subsequently exit from 
the resonance. As the exit also happens at a ran¬ 
dom value of the inclination cycle, the values of or¬ 
bital inclination with which bodies enter the cycle 
typically differ from those with which they exit. In 
this way the vi^ resonance changes the inclination 
of objects that cross it. Certainly, this mechanism 
would not work without the Yarkovsky effect. 


Interestingly, this process is very similar to the 
one observed inside another relatively week secu¬ 
lar resonance, namely the g + 2g^ — 3ge, which 
affects members o f the Koronis asteroid family 
( Bottke et al. [200lh . 


PLACE FIGURE 4 HERE. 


It is also worth mentioning that although pas¬ 
sage across the I'lc secular resonance could result 
in a very fast change in orbital inclination, the to¬ 
tal changes are limited and could not exceed the 
maximal variations. The amplitude of variations 
in orbital inclination is not the same for all ob¬ 
jects, but it is generally similar as the one shown 


in Fig. m 

Moreover, we have found that our test parti¬ 
cles typically spend 15-30 Myr inside the vic, be¬ 
fore being moved outside by the Yarkovsky effect. 
Thus, as a libration period is very long (about 
40 Myr), most of the particles spend less than one 
librating cycle inside this resonance (Fig. |3|). 

Finally, let us re-examine the role of zi for 
the members of the Hoffmeister family. The re¬ 
sults presented above clearly indicate that with¬ 
out Ceres in the model, the zi does not affect the 
family. Still, with Ceres included in the dynami¬ 
cal model zi may only affect a limited number of 
members. Our analysis has shown that about 5 
per cent of family members could reach this reso¬ 
nance after their inclination is pumped-up, at least 
a bit, by the vic resonance (see Fig. E]). The or¬ 
bital inclination of these objects is then addition¬ 
ally dispersed by zi, being the main reason for the 
slightly different distributions towards high and 
low inclinations. Moreover, although the analysis 
along this line is beyond the scope of the letter, we 
noticed that the zi resonance also slightly affects 
the orbital eccentricity of family members. 

Nevertheless, the role of the zi in the dynami¬ 
cal evolution of the Hoffmaister family members is 
minor compared to the role of the vic secular reso¬ 
nance. Thus, the evolution of the family is almost 
completely determined by the combined effect of 
the Vie resonance and the Yarkovsky effect, and 
would be practically the same even if the zi is not 
that close. 

3. Conclusions 

This is the first time a compelling evidence for 
orbital evolution of small bodies caused by a sec¬ 
ular resonance with an asteroid has been found. 
We prove that the post-impact transformation of 
the Hoffmeister asteroid family is a direct conse¬ 
quence of the nodal secular resonance with Ceres. 
This result has very important repercussions for 
our view of how Ceres-size bodies affect the dy¬ 
namics of nearby objects, and opens new possi¬ 
bilities to study such effects in the main asteroid 
belt and beyond. Examples include the dynamics 
of specific asteroid populations, the early phases of 
planetary formation, and extra-solar debris disks. 
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FIGURE CAPTIONS: 

Fig. 1 : The HofFmeister family in the space 
of proper orbital semi-major axis versus sine of 
proper orbital inclination. Note the strange shape 
of the family in this plane, in particular the large 
dispersion in the sin(ip) direction of the part lo¬ 
cated at semi-major axis less than about 2.78 AU 
(denoted by a shaded area in this figure). 

Fig. 2 : The evolution of the Hoffmeister fam¬ 
ily in the space of proper orbital elements. The 
three panels show the distribution of the test par¬ 
ticles after 5, 15, and 145 Myr of the evolution, 
from top to bottom, respectively. The orange dots 
represent the evolution of the particles within the 
dynamical model that includes the four giant plan¬ 
ets, from Jupiter to Neptune, and accounts for 
the Yarkovsky effect. The dark-green dots show 
the evolution when Ceres is added to the previ¬ 
ous model. The vertical dashed lines mark the lo¬ 
cations of mean motion resonances. The inclined 
blue lines denote the position of the zi secular res¬ 
onance; the solid line refers to the center, and the 
dashed line refers to the approximate lower bor¬ 
der of this resonance. Finally, the red dashed lines 
mark the approximate position of the vic = s — Sc 
secular resonance with Ceres. These plots clearly 
show very different dynamical evolutions of parti¬ 
cles, depending whether or not Ceres is considered 
as a perturbing body. While after 5 Myr the two 
distributions are still quite similar, after 15 Myr 
a remarkable difference is visible, with green dots 
evolving along the vic resonance. The last snap¬ 
shot corresponds to the distribution of the green 
dots after 145 Myr, very similar to the distribu¬ 
tion of the real Hoffmeister’s family members, in¬ 
dicating that this resonance is responsible for the 
strange shape of the family. 

Fig. 3 : Time evolution of the critical angle of 
the Vie = s — Sc secular resonance (top), and the 
mean inclination (bottom) for one of the test par¬ 
ticles. The correspondence of the time periods in 
which the critical angle is librating and the mean 
inclination is rapidly increasing (marked in both 
panels by the shaded area) clearly reveals the im¬ 
portance of the Vie secular resonance with Ceres 
in asteroid dynamics in this part of the main as¬ 
teroid belt. The solid black line shows the average 
of the mean inclination to better appreciate the 
evolution. 


Fig. 4: Variation of the mean orbital inclina¬ 
tion im vs. the critical angle — fie for the same 
test particle as shown in Fig. [3l. Note, however, 
that only time interval when this particle is inside 
the resonance is shown, i.e. from about 66 to 85 
Myr. During this time span the inclination of the 
particle undergoes cyclic variations with two dif¬ 
ferent periods. The first mode of these oscillations 
(red points) has a very short period of about 38 
kyr; if averaged out, it reveals the second, long 
period one (denoted by the black line), with a pe¬ 
riod of about 25 Myr. Note that the long period 
is actually associated to the libration of the crit¬ 
ical angle. A similar situation was also observed 
for secular resonances wit h the major planets (see 
e.g. iFroeschle et ah 1991 1. 
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